Structure and electronic properties of epitaxial graphene grown on SiC 
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We review progress in developing epitaxial graphene as a material for carbon electronics. In 
particular, improvements in epitaxial graphene growth, interface control and the understanding of 
multilayer epitaxial graphene's electronic properties are discussed. Although graphene grown on 
both polar faces of SiC is addressed, our discussions will focus on graphene grown on the (0001) 
C-face of SiC. The unique properties of C-face multilayer epitaxial graphene have become apparent. 
These films behave electronically like a stack of nearly independent graphene sheets rather than a 
thin Bernal-stacked graphite sample. The origin of multilayer graphene's electronic behavior is its 
unique highly-ordered stacking of non-Bernal rotated graphene planes. While these rotations do not 
significantly affect the inter-layer interactions, they do break the stacking symmetry of graphite. It 
is this broken symmetry that causes each sheet to behave like an isolated graphene plane. 



I. INTRODUCTION 

Although graphene was first isolated by chemical ex- 
foliation in 1961^ and later shown to grow epitaxially 
on SiC in 1975^ and 1998,1^ research on graphene's elec- 
tronic properties did not begin until 2001.HE1 As is so 
often promoted, the driving motivation for this research 
explosion is graphene's potential for carbon electronics. 
For this end game to be reachable, there are two im- 
portant properties of graphene that must be achieved. 



First, no matter how graphene is produced, it must have 
the band structure (and thus the transport properties) 
of an ideal isolated graphene sheet. Second any method 
of graphene production and device fabrication must be 
scalable up from a single nanometer scale prototype 
to macroscopic systems in which millions of graphene 
switches and their interconnects are integrated. Cur- 
rent research has focused on three ways of producing 
graphene: (i) epitaxial graphene grown on SiC,'^'^ (ii) 
graphene grown on metals,—'^ and (iii) exfoliated graphene 
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FIG. 1: 2D Brillouin zone of graphene near Ef showing the 
six Dirac cones at the A'-points. The cones are shown rotated 
through an angle (p relative to the SiC (2130) direction. 



separated from bulk graphite crystals and deposited on 
Si02 substratesl ^^ l ^^l The latter two methods both re- 
quire an end step: the transfer of graphene to a semi- 
conducting or insulating device platform. This transfer 
step is recognizable as the same problem that has plagued 
the development of carbon nanotube (CNT) electronics 
for the past two decades. In fact, the original impetus to 
develop epitaxial graphene electronics was to take advan- 
tage of CNT's unique electronic properties while finding 
a way to circumvent CNT's lack of scalability.'^^ The 
scalability of epitaxial graphene is well recognized. It has 
already been used to construct multi-FET array^^^ and 
has recently been used to build a 10,000 FET transistor 
array on a 3.5 x 4.5mm SiC substrate.'^ 

In addition to its scalability, epitaxial graphene 
beha ves electronically like an isolated graphene 
gj^gg^jmilHlII jts relevant electronic signature is the 
dispersion of the tt- and tt* -bands near the six if-points 
of the graphene hexagonal reciprocal unit cell.^^ The dis- 
persion is linear E{Ak) = hvp^k; where vp ^lO^m/sec 
is the Fermi velocity and Afc is the momentum relative 
to the -fC-points. The two-dimensional dispersion is 
isotropic and defines a cone (referred to as a Dirac cone) 
with an apex at the Dirac point, For undoped 

graphene the Fermi energy, Ep, coincides with Ejj 
so that the Fermi surface consists of six points [see 
Fig. [l]. This unique electronic structure is relevant 
for graphene based electronics for several reasons. For 
instance, the magnitude of vp means that electrons 
with energies significantly larger than thermal energies 
(~ leV) relative to Ed have wavelengths of the order of 
Airhvp /leY 2mn. Consequently quantum confinement 
becomes important in nanoscopic graphene structures 
with these dimensions and can lead to band gaps of the 
order of ^ leV.'^ In addition the symmetry change of the 
electron wave function as an electron moves from K to 
K' , is responsible for the ballistic transport properties 
of graphene .1221 

While most forms of graphene show many of the prop- 
erties of an isolated graphene sheet, only multilayer epi- 
taxial graphen e (MEG) g rown on the C-face of SiC ex- 
hibits them allEE^lMIl in addition. Landau level spec- 
troscopy from MEG films has demonstrated unprece- 
dented graphene properties including exceptionally high 



room temperature mobilities (> 250, OOOcm^/Vs), re- 
solved Landau levels in magnetic fields as low as 50 mT 
and remarkably low electron-phonon coupling up to room 
temperature.l^ 

Despite these impressive properties, it is often incor- 
rectly stated that epitaxial graphene is not true isolated 
graphene. The arguments imply that because the tt- and 
TT* bonds responsible for electron transport are perpen- 
dicular to the graphene plane, it is reasonable to con- 
clude that either placing graphene on any substrate or 
stacking graphene in multiple sheets will have a signif- 
icant impact on its band structure near the X-points. 
While this can be true, the literature oversimplifies the 
actual physics of the graphene-SiC substrate or even the 
graphene-graphene interaction. This thinking imposes 
unnecessary constraints on a number of potential re- 
search avenues and can actually hinder the pace of devel- 
oping graphene electronics. Therefore, before discussing 
the details of epitaxial graphene, it is worth taking some 
time to clarify these issues. 



A. Graphene-substrate interactions 

Graphene tends to interact more strongly with most 
materials than with itself. In fact this is the reason 
that graphene can be exfoliated onto a Si02 surface; 
the graphene adheres to the Si02 strong enough to al- 
low one or multiple layers to cleave from bulk graphite. 
The question is not if graphene bonds to a surface but 
how the bonding affects the electronic properties of the 
film. In the case of exfoliated graphene on Si02 the 
graphene substrate bond has not been extensively stud- 
ied. This is due in part to the problem of finding small 
fiakes with standard surface analysis probes. Even in 
the case where the band structure of exfoliated films 
has been directly probed by angle resolved photoemis- 
sion spectroscopy (ARPES), the film disorder severely re- 
duces the /c-resolution of the technique and makes quan- 
titative analysis problematic.'^ Transport measurements 
indirectly show the effects of the substrate: large spatial 
charge fluctuations and electron-doping levels as high as 
300meV.I^2l281 indeed many of these problems can be re- 
duced by an order of magnitude simply by removing the 
substrat e and suspending exfoliated graphene over etched 
channels .'22l2ni xhe change in transport and doping be- 
tween graphene with and without the Si02 substrate in 
part demonstrates the influence of the substrate interac- 
tion. 

However, not all substrate interactions are created 
equal. Understanding the differences in how graphene in- 
teracts with its environment, has important consequences 
in how to manipulate graphene's transport properties. 

For epitaxial graphene grown on either the SiC(OOOl) 
Si-terminated face (Si-face) or the SiC(OOOl) C- 
terminated face (C-face), sp'^ bon ding causes a strong 
graphene-substrate interaction.'^^Illn fact the interaction 
is so strong on the Si-face that the graphene interface 
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FIG. 2: Schematic models of graphene on the SiC Si-face. 

(a) Strongly bonded graphene zero-layer (ZL) or buffer layer. 

(b) graphene layer above the buffer, (c) H2 passivated buffer 
layer, (d) H2 passivated surface with two graphene layers. 
From Ref. 19 . 



the Si-fac^and 0.2eV on the C-fac^due to substrate 
charge transfer. It is not clear if the charging is from 
an intrinsic Schottky barrier, defect SiC danghng bond 
states at the interface or both.'22l It is clear, however, 
that these states can be removed by hydrogen passiva- 
tion. Riedl et a\^^ have been able to intercalate molecu- 
lar hydrogen between the Si-face SiC and the buffer layer, 
which simultaneously breaks the ZL graphene SiC sp^ 
bonds and saturates the remaining SiC dangling bonds 
[see Figs.[2]^c) and (d)]. These experiments are a beauti- 
ful example of how interface manipulation can eliminate 
the charge transfer to the first few graphene layers and 
restore the linear dispersion of the buffer layer in epitax- 
ial graphene [see Fig. |3jb)] . 
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FIG. 3: ARPES data near the A'— point of graphene grown on 
the Si-face of SiC. (a) for ZL graphene there are no observable 
bands, (b) Dirac cone from the ZL graphene after H2 passiva- 
tion, (c) Doped Dirac cone from a graphene layer above the 
ZL. (d) AB bilayer bands from ZL and the first graphene layer 
in (c) after H2 passivation of the interface. From Ref. . 



layer, sometimes called the zeroth layer (ZL) or buffer 
layer, becomes a wide gap semiconductor with none of 
the electronic character of isolated graphene [see Fig. [2]^ a) 
and Fig. |3][a)].'^ Nonetheless, the existence of this buffer 
layer has a major advantage; it largely decouples subse- 
quent graphene layers from the substrate. This has been 
demonstrated very nicely on Si-face epitaxial graphene 
where the linear bands of graphene are intact (although 
doped) in the layer above the buffer layer [see Fig. ^h) 
and Fig.[3]^c)].^3 

While the buffer layer preserves the linear bands, the 
graphene layers are still n-doped by as much as 0.44eV on 



B. Symmetry and graphene-graphene interactions 

The role of symmetry in determining the band struc- 
ture of graphene sheets can be demonstrated by compar- 
ing graphene with different stacking arrangements. In 
graphite the graphene sheets are rotated 60° relative to 
adjacent sheets in the stack (Bernal stacking) .^^ Bernal 
stacking causes the two atoms per graphene cell to be in- 
equivalent (the two atoms are labeled 'A' and '5'). The 
atoms are bonded to an atom in the graphene sheets 
above and below through orbitals. The '5' atoms lie in 
the center of graphene hexagons in sheets above and be- 
low and are therefore not bonded to other carbon atoms. 
The stacking induced in-equivalence of the and '5' 
atoms in graphite is often referred to as AB stacking. 

Si-face graphene is Bernal stacked like graphite. '^'^ This 
has been elegantly demonstrated by comparing Si-face 
graphene that has been grown with both a ZL and an ad- 
ditional graphene layer before and after H2 passivation.l^ 
In the as-grown Si-face graphene, only the top layer is iso- 
lated and exhibits the band structure of a single graphene 
sheet [see Figs. [Sf^c)]. When the H2 is intercalated into 
the interface, the ZL becomes isolated from the substrate 
and the two graphene layers become a new electronic sys- 
tem. Because the ZL is rotated 60° relative to the top 
layer (Bernal stacking), the doped Dirac cone of the sin- 
gle layer becomes the split bilayer bands structure of a 
bilayer pair [see Fig. [sj^d)]. 

While it is tempting to assume that the band struc- 
ture change is due to the weak bonding of the tt orbitals 
between planes, this is incorrect. The change is instead 
due to the symmetry change of the bilayer system. To 
see this we can examine the band structure of bi-layer 
graphene when interplanar interactions exist but the AB 
stacking symmetry is destroyed. The simplest example is 
AA stacking.El In AA stacking (0° rotation) the number 
of bonds/area is double that of Bernal stacking (one for 
every atom). Despite this increase, the dispersion in A A 
stacking reverts back to the linear bands of graphene.!^ 
This effect is more general and applies to any rotation 
other than 60°. 

Although non-60° rotations lead to small increases 



FIG. 4: A model of a 7x7 commensurate structure produced 
by two graphene sheets with a relative rotation of 21.79°. 
The top sheet is represented by the ball and stick model. 
The lower sheet is represented by open and grey circles. In 
the lower layer A atoms (open circles) and B atoms (grey 
circles) are directly below and bonded to atoms in the upper 
layer (shaded red or blue) with equal probability. Although 
not easily visable, shaded and open atoms in the layer below 
occupy positions in the upper lattice with equal probability. 



FIG. 5: Calculated band structure for three forms of 
graphene. (i) isolated graphene sheet (dots), (ii) ^45 graphene 
bi-layer (dashed line) and (iii) a {y/l3x \/13)gR-46.10° struc- 
ture formed by two stacked graphene sheets with a relative 
rotation angle of 32.204° (solid line). Note that subscript "G" 
refers to graphene lattice vectors. Inset shows the details of 
band structure at the Tf-point. From Ref. |38| . 



in the bonding energy per atom (of the order of a 
few meV/atom),"^^ these small energetics concerns are 
overshadowed by a more important symmetry change 
that leads to dramatic changes in the electronic band 
structure of rotated graphene sheets. Here rotation 
means that the relative angle between adjacent sheets 
is some value other than 60°. To demonstrated this, 
we show a ball model of two graphene sheets rotated by 
21.79° in Fig.l4](a similar structure can also be formed by 
rotating 16.43°). For the commensurate structure shown 
in Fig. |4j the in-equivalency of the and 'B' atoms, as- 
sociated with Bernal stacking, is not longer true. There 
are as many atoms in sites below atoms in the up- 
per plane (bonding sites) as there are '5' atoms in the 
same position. Actually, all 'j4' atoms in the layer below 
occupy positions in the upper graphene lattice equally 
likely as 'B' atoms. In other words there is no symmetry 
breaking in the sheet so and 'B' atoms are essentially 
equivalent in terms of bonding and structure. This result 
is true for any rotation angle other than 60°. The effect 
of this symmetry is dramatic. Both planes become elec- 
tronically equivalent to an isolated graphene sheet. This 
has been demonstrated theoretically a number of ways 
for both small and large relative rotation angles .'33321 
Figure [5] compares the calculated band structure for two 
graphene sheets rotated by 32.204°. The electronic bands 
are obviously different from the AB stacked pair and in- 
distinguishable near the Dirac point from the bands of a 
single graphene layer. 

In C-face graphene the majority stacking is non- 
graphitic so that films as thick as 60 graphene layers 
still behave electronically like a stack of isolated graphene 
sheets. f^^Hi! In Scc.|ll]we wiU show that AB stacking is the 
exception rather than the rule in C-face films and that 
the band structure of these multilayer films are equivalent 



to isolated graphene. 

II. STRUCTURE AND STACKING 

Graphene grown on the C-face of SiC is significantly 
different than Si-face graphene, both in how it grows and 
its structural order While the slow growth rates of Si- 
face graphene makes it relatively easily to grow 1-3 layers 
films, the high growth rates in C-face graphene makes 
thin film growth much more difficult. In addition. Si- 
face epitaxial graphene grows Bernal stacked (i.e, with 
a 60° relative rotation between adjacent planes), while 
C-face graphene grows in an ordered set of relative rota- 
tional angles. In this section we will focus on the detailed 
morphology and electronic band structure of C-face mul- 
tilayer epitaxial graphene films. In Sec. |III| we will show 
that the rotational stacking in C-face graphene is highly 
ordered and is not disordered as is often asserted in the 
literature. 



A. Structure of C-face graphene 

Graphene grown on the Si-face is always oriented 
30° relative to the SiC (2130) direction and leads to a 
(6\/3 X 6\/3)R30° reconstruction in low energy electron 
diffraction (LEED) images [see Fig. |6|a)].^ Graphene 
grown on the C-face is also rotated 30° relative to the 
(2130) direction but LEED and surface x-ray diffraction 
(SXRD) indicates that some planes are rotated within a 
small angularregion around the SiC (2130) [see Fig.rolb) 
and Fig.g.l^Sl! 

The additional rotation angles are seen in both Ultra 
High Vacuum (UHV) and furnace grown graphene. How- 
ever, it is important to recognize that there are significant 
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FIG. 6: LEED patterns from (a) Si- and (b) C-face graphene 
films grown on SiC. Both images show SiC diffraction spots. 
In (a) the graphene is rotated 30° relative to SiC and the 
(6%/3x6-\/3)R30° reconstruction spots are visible. In (b) a set 
of six diffuse arcs rotated 0° from SiC are also visible. 




FIG. 7: Schematic LEED pattern from a C-face graphene 
film. Both graphene spots rotated 30° from the SiC (2130) 
direction and diffuse graphene arcs centered around ^=0° are 
shown. A hexagonal graphene Brillouin Zone (BZ) rotated by 
(j) is shown. Note that the FK direction of the graphene BZ is 
rotated 30° relative to the graphene reciprocal lattice vector 
ah- 



morphological differences in the graphene sheets grown 
in UHV and at ambiant pressure furnace. In Fig. [8] we 
compare low energy electron microscopy (LEEM) micro 
LEED images of UHV and furnace grown graphene. The 
image for a UHV film shows both extra rotated graphene 
spots and diffuse graphene arcs. This is evidence of 
a large number of rotation angles within the ^ 1/im 
beam diameter, consistent with STM studies on UHV 
samples.!^ Furnace grown samples, on the other hand, 
show only a few different graphene rotations in the same 
area. This emphasizes that while both methods produce 
rotat ed she ets, the size of the sheets in UHV films is very 
small.l^i'^ In the furnace grown samples, the size of the 
graphene sheets is much larger. In fact, grain boundaries 
have yet to be seen in furnace grown graphene, suggesting 
that continuous graphene sheets can span macroscopic 
dimensions. 

It is very important to realize that the C-face rotated 
diffraction patterns, like the one shown in Fig. [6][b), 
are not due to randomly rotated domains of graphite 
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FIG. 8: LEEM micro LEED images of (a) a UHV grown 3- 
to 4-layer C-face graphene film. Only one quadrant of the 
pattern is shown but the two primary graphene reflections 
are given for reference. Note diffuse arcs and extra graphene 
spots, (b) a furnace grown 11-layer graphene film on the C- 
face of SiC. Only two rotated graphene planes are visible. The 
image is slightly distorted from aberrations in the LEEM. 

crystals like those in HOPG graphite. Instead, the ro- 
tated sheet are interleaved in the multilayer graphene 
stack .ti^ilMl This statement is the result of a wide va- 
riety of experiments including transport measurements, 
electronic structure and structural studies. These will 
be discussed in detail below and in Sec. IIIIl Because of 
C-face graphene's unique electronic properties and scal- 
ability, we will focus on furnace grown C-face graphene 
in the rest of this review. 

To begin to understand the rotational structure and 
its implication to C-face graphene's electronic proper- 
ties, we begin by correlating SXRD and ARPES data. 
To make this easier Fig. [7] shows a schematic LEED pat- 
tern from a C-face graphene film like the LEED pattern 
in Fig. [6][b) . For reference a graphene sheet (with a re- 
ciprocal lattice vector Oq) rotated by an angle relative 
to the SiC (2130) produces diffraction intensity on the 
graphene arc at the same angle 0. The corresponding 
graphene Brillouin zone (BZ) is rotated by </> so that the 
YK direction vector is rotated by an angle a = + 30° 
relative to the SiC (2130) direction [see Fig. [?]. 

Figure [9] shows the angular distribution of graphene 
planes near = and 30° for two different samples. 
There are three things to point out in the figure. First, 
as seen in the LEED, the angular distribution near = 
is much broader than the distribution near = 30 (note 
the expanded 0-scale in Fig. |9](b)). The reason for this 
difference becomes obvious by considering the possible 
graphene structures that are nearly commensurate with 
SiC (0001) when a rotated graphene plane is placed on 
the surface.^ We have marked the angles for all (L x V) 
and (L-\/3x-£a/3) R30° graphene-SiC near commensurate 
structures in Fig. [9] For comparison, the height of the 
drop lines for the markers has been scaled by the inverse 
lattice mismatch between the graphene super cell and 
the SiC substrate: 1/AL. This way of plotting accentu- 
ates structures with both small strain, e, and small cell 
sizes (AL = eL). In addition to graphene being nearly 
commensurate with the SiC, two graphene sheets can be 
rotated relative to each othei"^* to form a commensurate 
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FIG. 9: X-ray angular distribution of graphene rotations near 
(a) (f) = 30 and (b) (j) = for two different C-face samples. 
Sample #f (□) is a 30-layer film and (o) sample #2 is a 6- 
layer film. Horizontal dashed line shows the offset for sample 
#f data. Drop lines show the position of graphene-SiC near 
commensurate structures for (red A) LxL and (blue A) L\^x 
I/\/3 cells. The height of a drop line is 1/AL, where AL is 
the mismatch between the graphene structure and the SiC 
lattice. (Grey o) show the position of rotated graphene sheets 
that are commensurate with a graphene sheet rotated (p = 30° 
relative to SiC (drop line height is proportional to f/C for a 
^/CxVCRe unit cell). 



(■v/C* X \^)qRO graphene-graphene super cel\^ Figure 
^h) marks the angular position of all graphene sheets 
that form a super cell with a graphene sheets rotated by 
(j) = 30° with respect to the SiC. The height of the drop 
lines for these graphene-graphene commensurate angles 
in Fig.jojhas been drawn proportional to 1/C to highlight 
small unit cell structures. 

It is clear that the density of graphene-SiC near com- 
mensurate rotations per unit of arc is higher for the 0° 
rotations than for the 30° rotations. In other words there 
are simply more possible commensurate structures ori- 
ented around 0°, with small energy differences."^^ The 
higher entropy associated with the distribution of com- 
mensurate angles near = 0° would explain why the 
SXRD angular distribution is broader in the (f> = az- 
imuth. 

The second thing to note in Fig. [9] is that the exact dis- 
tribution or rotation angles is sample-dependent. There 
are a number of possible reasons for this that can be 
related to the small energy differences for different ro- 
tational angles. It is known on at least the Si-face that 
graphene grows out from substrate step edges.*^ This 
means that a specific graphene orientation could be influ- 



enced by slight rotations as the graphene grows from the 
SiC step. Since the substrate step direction depends on 
sample miscut and polishing, the distribution of graphene 
orientations would be sample dependent. The pleats in 
graphene that form when graphene is cooled from the 
growth temperatur^^ can be another source of small ro- 
tations that do not introduce defects in the hexagonal 
lattice of an otherwise continuous film. Scanning tunnel- 
ing microscopy (STM) has shown that rotational changes 
do occur at these topological boundaries.l^Sl Regardless 
of how the distribution is initially set, the rapid quench 
when the graphene is cooled from the growth tempera- 
ture would freeze in the distribution. Whether or not 
these rotations can be annealed out remains to be seen. 

Finally, a detailed analysis of the SXRD angular dis- 
tribution shows that the integrated area of the (f> curves 
around = 30° is nearly the same as the area around 
0° (i.e J hodcj}/ J lodct) ^ 1± 0.2). In other words it is 
nearly as likely to find a graphene plane rotated 30° rel- 
ative to SiC as it is to have one rotated 0°. This implies 
some order to the rotational stacking sequence and will 
be discussed in more detail in Sec. IIIII 

Why the rotated graphene planes around 0° form on 
the C-face graphene and are not produced during Si- 
face growth is not understood. It has been argued that 
the coupling between SiC and the first graphene layer 
is stronger on the Si-face compared to the C-face.'^^'^ 
The existence of a stronger Si-face interaction is used to 
explain the observed difference in angular distribution 
on the two SiC surfaces. The argument suggests that 
the stronger Si-face interaction forces the graphene to be 
aligned ±30° from the SiC, while the weaker C-face in- 
teraction allows the graphene to orient itself in multiple 
rotation directions. Other groups, on the other hand, 
have argued for a stronger C-face interaction compared 
to the Si-face based on both graphene-SiC bond lengths 
and inverse photoemission data.l^SlIlUlIl 

These arguments may be missing the point, because 
stronger or weaker interactions do not explain the differ- 
ent stacking order on the two surfaces. In particular, a 
strong Si-face interaction does not explain why graphene 
grown on the Si-face is AB stacked. Regardless of the in- 
teraction strength, the second layer must "know" that it 
should be rotated 60° relative to the first. It could be ar- 
gued that the AB stacking pattern is formed by high tem- 
perature annealing of rotated sheets to the lower energy 
Bernal state. However, a similar ordering should, by the 
same logic, occur in C-face films that are grown at higher 
temperatures. Since this does not happen, we are forced 
to propose something about the interface template that 
forces the AB stacking order in Si-face graphene. Simi- 
larly, a strong C-face interface interaction could explain 
why the graphene rotations angles in general match low 
strain and small unit cells near commensurate graphene- 
SiC structures. However, it does not explain why C-face 
graphene grows with a more ordered rotational stacking 
sequence where planes are on average 30° apart (as we'll 
see in Sec. HI ). As in the Si-face case, a symmetry change 
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FIG. 10: ARPES scans taken at the if-point radius {k^ = 
1.704A) for a 10-layer grpahene film on the C-face of SiC. The 
photon energy is 36eV. The scans are taken at two different 
emission directions: (a) along the SiC (2130) (a = 0°) and 
(b) (1010) (q = 30° 
Fig. [3 



must occur at the interface that cause the 30° rotations 
to develop. This view makes the important question not 
which surface has the stronger interaction, but what is it 
about the different interface structures that imposes dif- 
ferent rotational stacking sequences for graphene grown 
on the two SiC surfaces. 



B. Band Structure of C-face graphene 



As discussed in Sec. IB the rotational stacking is ex- 
pected to preserve the symmetry of an isolated graphene 
sheet. Experimentally this means that the band struc- 
ture of multilayer epitaxial C-face graphene should con- 
sist of a large number of multiply-rotated graphene Bril- 
louin zones with the same distribution of rotation angles 
as seen in the SXRD data. Figure [TO] shows ARPES 
scans at the ii'-point radius {k^ — 1.704 A) for elec- 
trons emitted along the SiC (213Q and (1010) directions. 
Both scans show multiple linear dispersing Dirac cones 
but the distribution of these cones in the (1010) direc- 
tion is bimodal and peaked at Aky ~ ±0.2A ^ which, 
for small rotation angles, correspond to cones rotated 
a = 30° ± Aa«30° ±tan-i(Afcy/fcrif) = 30°±6.7°. Be- 
cause the graphene TK direction in ARPES is rotated 
30° from the graphene reciprocal space direction, Gq [see 
Fig. [t], the split distribution of graphene rotation an- 
gles observed in LEED around cf) — would produce a 
split distribution of if -points rotated around a = 30, i.e. 
the Dirac cone distribution would be centered around 
the (1010) direction and peaked at ^6 — 7° as observed. 
Likewise the narrow distribution of graphene rotations 
around (j) — ^30 would produce a narrow distribution of 
Dirac cones along the SiC (213(| (a — O) direction as seen 
in Fig. [To] 

To emphasize the correlation between graphene rota- 
tion angle </> and the TK rotation direction a, we have 
marked the discrete rotation angles, a, of the ARPES 
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FIG. 11: A comparison of the graphene SXRD angular dis- 
tribution in (j) (o) and the ARPES Dirac cone histograms in a 
(Solid lines), (a) SXRD graphene distribution taken around 
= and corresponding ARPES histogram of cones dis- 
tributed around a = 30. (b) SXRD distribution taken around 
(ji = 30 and corresponding ARPES histogram of cones dis- 
tributed around a = 60. Insert in (a) shows a magnified view 
of a single rotation angle. 



Dirac cones from the top three graphene layers against 
the bulk angular distribution of graphene rotations, (j), 
measured by SXRD in Fig. 



11 It is clear that the 



ARPES cone positions correlate well with the data for 
graphene rotations in both the = 30° and 0° direc- 
tions. Two experimental differences between SXRD and 
ARPES should be noted. First the ARPES is only mea- 
suring cones in the upper 1-4 layers, while SXRD is mea- 
suring graphene planes throughout the entire film. Sec- 
ond the SXRD beam size is 3mm while the ARPES beam 
size is 40/j,m; this is why ARPES data shows a small num- 
ber of discrete rotated cones and SXRD shows a more 
continuous distribution averaged over a large beam foot- 
print. Also note that the angular width of each discrete 
rotation is very narrow; a detailed scan of one such angle 
is shown in the insert of Fig. 11 a). Its width is 0.045°, 



corresponding to an x-ray rotational coherence distance 
of ^ 1/im. This simply confirms that the x-rays can only 
measure graphene lengths up to the distance between SiC 
steps 1/im for these samples) after which x-ray diffrac- 
tion looses coherence as the graphene flows over the SiC 
substrate steps.'^ 
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FIG. 12: Three rotational stacking models, (a) Bernal stack- 
ing sequence with 50% randomly placed rotational planes 
(dashed lines). Arrows mark isolated A or B type graphene 
sheets, (b) Quasi-ordered stacking with the same concentra- 
tion as (a) except the rotational plane sequence is nearly or- 
dered, (c) Interrupted stacking sequence where the rotational 
planes cause a random phase boundary (marked by brackets) 
in the Bernal stacking sequence. Upper grey pattern in (c) 
represents the unperturbed sequence in (a) to show the dif- 
ference between random stacking and the interrupted model. 



The discrete rotated set of Dirac cones clearly demon- 
strate that the rotated graphene planes have become 
electronically identical to isolated graphene sheets with 
nearly perfect linear dispersion at the X-point. Be- 
cause the C-face films are thick, charge transfer from the 
SiC interface to the upper graphene layers is insignifi- 
cant. However, doping is observed in the graphene sheets 
at the graphene vacuum interface. Doping is predomi- 
nantly p-type but can vary to n-type depending on the 
sample.^ The reason for this sample-dependent doping 
is not known although weakly bound adsorbates could 
be the source. The doping ranges from ~ 33meV p- 
doped on some samples to n-doped as low as -14meV on 
others. This gives a charge density that ranges between 
~ 10^^ — 10^°cm~^, comparable to IR measurements from 
similar films (5 x lO^cm"^).'^^ 



III. MEASUREMENTS OF THE STACKING 
ORDER 



are randomly introduced into an otherwise AB stacked 
film with a probability 7. Figure |12[ a) shows a ran- 
dom sequence of rotational planes. Note that an A plane 
must be adjacent to at least one B plane to locally have 
the band structure of an AB pair [marked by brackets 
in Fig. |12[ a)]. A or B planes surrounded by rotational 
planes, [marked by arrows in Fig. |12[ a)], will still have 
the band structure of an isolated graphene sheet. For the 
random model the probability that any plane is part of 
an AB pair is determined by the two ways of retaining 
an AB pair. This is either an ''..ABR." or ''..ABA." 
sequence. Pab is then the sum of the probabilities that 
these two sequences exist: 

Pab = 27(1 - 7)' + (1 - 7)' = (1 - 7)'(1 + 7) (1) 

In the random model there is always a chance that an A 
or B plane is surrounded by two rotated non-^B planes, 
even for very small 7's so that Pab is always less than 
1-7- 

However, once the rotational plane stacking becomes 
ordered, Eq. [T]is no longer valid and it will overestimate 
Pab- In fact the overestimation can be very consider- 
able. Consider the random model case with 7 = 0.5. 
Then on average there is a rotated plane every 2 graphene 
planes and Pab — 0.375. However, if the distribution of 
rotational planes was completely ordered (i.e a rotation 
exactly every other plane as in Fig. [l2|b)), there would 
be no AB pairs, since each A or B plane would be sur- 
rounded by a rotated plane. 

In this section we discuss both x-ray diffraction and 
ARPES experiments that directly address the question 
of how the rotational planes in the MEG films are dis- 
tributed. Early x-ray reflectivity measurements esti- 
mated the rotational plane probability to be approxi- 
mately 7 ^ 0.4 based on an average expansion of the 
interplanar graphene lattice constant caused by tt* bond 
interference in a random stacking model.!^ This value of 
7 would correspond to Pab '^0.50 in the random model. 
As we'll show in this section, this is a gross overestimate 
of the number of AB planes because the rotational stack- 
ing is more ordered than initially assumed. In fact the 
ordering is high enough that the AB planes, rather than 
the rotational planes, can be considered as faults in the 
stack. 



So far we have shown that rotated graphene planes 
exist that have the band structure of an isolated graphene 
plane at the i^-point. We have yet to discuss how these 
rotations planes are ordered. As we'll demonstrate, C- 
face graphene has an extremely high rotational ordering 
that makes it a new carbon allotrope. 

If one wishes to understand the transport in MEG 
films, it is the stacking order, and not just the concentra- 
tion of rotational planes, that controls the fraction, Pab, 
of planes in the film that are part of an AB stack. To see 
why this is so, we can compare two simple stacking mod- 
els. One model is to assume that rotational planes, R, 



A. X-ray analysis 

In order to demonstrate how SXRD data can be used 
to shed light on the stacking order in C-face graphene 
films, we begin by writing the general scattered x-ray 
amplitude from an N-layer graphene film as a function 
of the momentum transfer vector q; q = ko — (where 
ki and ko are the incoming and scattered x-ray wave 
vectors, respectively). 



N-l 



A(q||,L) (X e 



'(o-SiCTi/cc)^ 



Lck/ca 



(2) 



k=0 
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The momentum transfer vector normal to the surface, 
g^, has been written in terms of a variable L defined as; 

— ttL/cg where cq — 3.35A is the graphite inter-layer 
spacing. Ck is the vertical position of the k^^ graphene 
sheet. Fk{q\\) is the form factor for the fc*^ graphene 
plane and depends on the rotation angle of the sheet. The 
gaussian term in Eq. [2] accounts for the graphene RMS 
roughness, aa, that res ults f rom the graphene draping 
over SiC substrate steps.'^^'^ 

For the (lOL) crystal truncation rod (CTR),'22jthe mo- 
mentum transfer vector parallel to the surface is \q^\ \ = 
IoqI =27r/aG(-y3/2). For the following discussions we as- 
sume that aQ is pointing along the SiC (1010) direction 
so that = ±30. i^fe(aQ) can then take three possible 
values: it is equal to Fa or Fb, the form factors of ei- 
ther an A or B graphene sheet corresponding to sheets 
rotated </> = 30° and —30°, respectively, or Fk{aQ) — 
corresponding to sheets rotated by 7^ ±30°. 

While a general analytic solution to Eq. [2] is not pos- 
sible for the (lOL) CTR, it can be solved for the case of 
Bernal stacking with random rotations like the model in 
Fig. 12 ^a). The probability of a non-Bernal rotation is 
defined as 7. Note that there are three such types of pos- 
sible rotations; an A plane following another A plane, a B 
plane following another B plane, and any plane that that 
is not an A OT B plane. The first two contribute intensity 
to the (lOL) CTR and the latter does not. We define the 
probability that a non-Bernal plane is an A or B fault as 
/3. Using these definitions, the average intensity I — AA* 
from Eq. |2] becomes; 
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FIG. 13: The SXRD graphene (lOL) CTR from a 30-layer C- 
face graphene film. The data (•) has been normalized by the 
surface roughness term in Eq. 3a using (j = 0.54A. The best fit 
random model X-ray reflectivity (line) using the 7 = 0.425, 
/3 = 0.26, iV = 18 is shown. 



even is; 



I {L) / cxp[~2{asic'^L/ cg)^] — B for L odd and 



/p(i = odd) Jl-7)2+7/?[7/?-2(l-7)] 



IpiL^even) (1 - 7)2 + 7/3[7/? + 2(1 - 7)] 



(4) 



^(qii 



L) oc e 



-2((TsiC7ri/cG) 



B 



sin^{TTLN/2) 
sin^(7ri) 



For random rotational planes in the AB stack, this ratio 
is 3 for any value of 7 as long as /3 = 0. In other words, the 
random model requires AA or BB fault pairs to change 
the /p(i = odd)//p(L = even) ratio. 

With these properties of the random model in mind, 
we can compare the experimental C-face graphene (lOL) 
CTR with Eq. [3a] Figure [13] shows the experimental C- 
face graphene (10£) intensity. The data has been divided 
by exjp[—2{asicT^L/cG)^] using crsiC =0.54A to correct for 
the SiC step roughness. The obvious difference between 
the data and the fit is the significant background from the 
random model. The magnitude of the background results 
from requiring the fit to match the experimental ratio of 
C(7, (3) = [(1 - 7)^ + 7^/3^ - (1 - 7)/?7]. (3c) ^^'^ odd and even Bragg peaks /p (odd) //p (even) = 1.4. 



x{2C - [C - 3(1 - 7)7/3] cos(7rL) 

+V3 [7^/32 - (1 - 7)2] sin(7rL)}) , (3a) 



i?(7,/3)=7^[(l-7)+/3(l-27)], 



(3b) 



The intensity consists of a background term B{'j, /?) and 
a set of Bragg peaks at L =integer defined by the sine 
function in Eq. [3a] 

There are three features of this model that should be 
pointed out. First, there is a finite background that 
scales with N and goes to zero when 7 = 0. Second, 
note that the full width at half maximum (FWHM), 
AL, of the Bragg peaks is independent of 7 or a and 
only depends on the number of graphene layers; AL ~ 
2/N. The background and the 7 independent diffrac- 
tion widths are common results for diffraction from ran- 
dom defects.'^ Finally, the random model predicts a spe- 
cific relationship between the intensity of the L = even 
and L = o dd B ragg points. If we normalize the inten- 
sity in Eq. 3a by exp[— 2(crgic7rL/cG)2], then the ratio 



of the normalized peak intensity above the background 



We note that there are a range of 7 and /3 that give simi- 
lar fits to the ratio with the same high background. The 
fit in Fig. JTS] was done using 7 = 0.425, and /3 = 0.26. 
The value of 7 was chosen to be consistent with pre- 
vious x-ray specular reflectivity estimates based on the 
interplanar graphene expansion for a random model.^ 
Regardless of the exact parameters in the fit, the small 
experimental background is the first indication that the 
rotational stacking is not random. 

An even more significant problem with the random 
model is the number of graphene layers required to fit the 
data. The best fit to the data in Fig. JT3] that reproduces 
the Bragg peak widths is N — 18. However, ellipsometry 
measurements on the same film estimate N to be nearly 
twice as large (30 ± 3). Likewise, x-ray specular reflec- 
tivity (gy = in Eq. pi, which is sensitive to the total 
number of graphene planes regardless of the rotational 
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FIG. 14: Comparison of diflterent model fits to the data. The 
dashed (red) line is a fit using random rotations {w = 0) and 
interruption of the AB stacking (a>0) as in Fig. 12 c). This 
best fit uses A'^ = 33, 7 = 0.545 and a switching probability of 
a = 0.1. Solid (black line is the best fit using a quasi ordered 
model with switching. The fit uses A'^ = 30, 7 = 0.55, «; = 0.05 
and a switching probability of q = 0.09. 



distribution, measures a film thickness of = 33 ± 2. 
This inability to fit the Bragg widths reflects the random 
models explicit lack of correlations in the rotational plane 
stacking. 

In order to go beyond the random model and begin 
to understand the stacking in multilayer C-face films, 
requires the ability to compare the experimental data 
with other more complicated models where analytic so- 
lutions are not possible. To do this we can generalize 
the CTR intensity calculations to a Quasi-ordered stack- 
ing (QOS) model. In this model the stacking sequence 
is generated by assuming two rotational planes are on 
average separated by I/7 (the same as in the random 
model). Disorder is introduced by allowing the rota- 
tional planes to occupy sites that are An from the average 
according to gaussian weighted probability distribution; 
oc exp[— (An— l/7)^/2w^]. The width of the gaussian (w) 
sets the degree of order: w = for perfect ordering and 
w = 00 for random stacking. To be completely general 
the model also allows for another type of stacking dis- 
order in addition to the rotational stacking. We assume 
that the AB sequence can be disrupted after a rotational 
plane is introduced. This is done by allowing a switch 
in the AB sequence from ..ABRBA... to ..ABRab... as 
illustrated in Fig. IT2jc) . The switch occurs with a ran- 
dom probability, a (a = means no switch occurs). Once 
a particular stacking sequence is set, Eq. [2]is calculated 
and the CTR intensity is calculated. The final calculated 
intensities are formed from an ensemble average of 1000 
randomly generated distributions. 

The interrupted stacking sequence is a logical exten- 
sion of a purely random model. A switch in the AB or- 
der will introduce phase boundaries in the stack that will 
broaden the Bragg peaks. A fit to the graphene (lOi) 



data using the interrupted model is shown in Fig.[T4j The 
effect of switching does indeed broaden the L=odd Bragg 
peak width so that even a = 33 layer films (consistent 
with the value determined by SXRD specular reflectiv- 
ity) can match the experimental width for a reasonable 
value of 7. Also note that unlike a purely random model, 
the ratio of odd and even Bragg peak intensities fit the 
experimental data without the need for any AA pairs (i.e. 
/3 = 0). Beyond these points the random model, even with 
switching, still predicts a significant background not seen 
in the data. In addition the width of the L=odd Bragg 
peak widths in this model remains a function of N only, 
just like the random model. It is clear that a more fun- 
damental change to the stacking order is needed to fit 
the data. 

Of the three models described in Fig 
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the quasi- 
ordered stacking model with a small amount of AB 
switching gives the best fit to the experimental data. As 



shown in Fig. 14 the model reproduces two features of 
the data that the random model cannot (even with AB 
switching); a low background and significant broadening 
of both the odd and even L Bragg peaks. 

In order to fit the low experimental background the 
order parameter was set to nearly zero w = 0.05, with 
the average spacing between rotational planes set to 
1/7 = 1.83. In order to fit the i = 1 Bragg rod, N 
was set equal to 30. This is only slightly smaller than 
the thickness derived from the specular reflectivity data. 
While the L = even Bragg peaks have broadened signif- 
icantly, the calculated rods are still narrower than the 
data. This is because some degree of randomness is left 
in the model. The sharp gaussian distribution set by 
ui = 0.05 and the non- integer value of I/7 requires that 
rotated planes landing at non-integer plane separations 
be randomly shifted to integer values (weighted to keep 
the rotational density equal to 7). 

The most important results of this analysis is the low 
value of the order parameter. It implies that the stack- 
ing sequence in multilayer epitaxial graphene is driven by 
predictive rather than stochastic processes. Presumably 
the SiC interface structure changes roughly every three 
SiC bilayers (the number of bilayers required to release 
the carbon necessary to form a single graphene layer). El 
This structural change forces a rotation of a forming 
graphene sheet to be 30° from the previously completed 
sheet. 

The low order parameter also predicts a low fraction of 
AB stacked planes. Figure [Ts] shows how the AB fraction 
depends on w. For 7's 0.5 and above, rotational stacking 
order causes Pab to fall to zero as w 0. The best fit 
to the graphene (lOL) rod data using the quasi ordered 
model with 7 = 0.545 leads to a value of Pab — 19.8%. 
This low value of Pab shows that rather than thinking 
of the rotated planes as faults in an otherwise Bernal 
stacked film, MEG films should be viewed as rotationally 
ordered non-Bernal stacked films where the AB planes 
are themselves the faults. What is surprising is that, 
despite their small concentration, the AB planes in the 
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FIG. 15: A plot of Pab versus the order parameter w for the 
quasi ordered model. The solid lines were calculated for a 
10-layer film with a = (no AA stacking) and /3 = (no AB 
switching). The dashed line shows the effect of A A stacking 
on Pab for the case of 7 = 0.6. 

film are placed as if they were in a Bernal stacked film, 
i.e. a is small. This again suggests a periodic graphene 
formation mechanism at the SiC interface. 



B. ARPES analysis 
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FIG. 16: Comparison of the ARPES band structure near the 
if-point for a single graphene sheet and an AB bilayer. The 
calculated tight binding dispersion from both (a) an isolated 
graphene sheet and (b) the band splitting from an AB stacked 
bilayer pair are shown, (c) and (d) two experimentally mea- 
sured bands for multilayered C-face graphene rotated = 30° 
from the SiC (2130). (c) shows only linear graphene bands 
while (d) shows both linear bands and a split band associated 
with an AB stacked bilayer pair. The photon energy is 36eV 
and the photon beam size is 40/im. 



ARPES data can also be used to corroborate the high 
degree of order inferred from the SXRD fits presented in 
the last section by estimating the relative number of AB 
planes in a C-face film. To understand how this is ac- 
complished, we first review how the AB planes are iden- 
tified and then discuss the experimental limitations of 
such measurements. 

The electronic signature of AB stacking is the split- 
ting of the linear bands at the i^T-point.'^^l This is shown 
in Figs. 16 |a) and (b). When two graphene sheets are 
Bernal stacked, the single Dirac cone splits into two 
parabolic sub-bands. The lower band is shifted nearly 
0.5eV to lower binding energies. It is important to un- 
derstand, as explained by Orlita et al.,'23 that the rela- 
tive strength of the two bilayer bands is modulated by 
the perpendicular momentum transfer, fcj^, defined as 
A:_L = (2m//i^)Vi?kin + Vq. Here -Ekin is the measured 
electron kinetic energy and Vq is the graphene inner po- 
tential, Vb = 16.5eVp2l Ohta et al.,'22l have shown that 

both sub-bands modulate with a period of fc_L'^2.0A , 

which is not 27r/cG = 1.86A as expected. The relative 
strength of the upper sub-band (compared to the lower 

band) goes to nearly zero every lA ^ . 

Figures [iBjjc) and (d) shows experimental cuts perpen- 
dicular to the YK direction (aligned with the SiC (1010)) 
to show the band structure from graphene planes rotated 
by (/) = — 30°. Both directions show the typical series of 
linear dispersing Dirac cones. In addition to the linear 
Dirac cones. Fig. [l6|c) also shows the lower parabolic 
sub-band that corresponds to an AB stacked graphene 
bilayer. Because there are so many cones closely spaced 
in ky in the C-face films, distinguishing between linear 



and parabolic bands can be difficult. Because the lower 
sub-band is much lower in energy, it is much more easily 
distinguished from linear cones than the upper sub-band. 
The problem of identifying AB bilayers can therefore be 
minimized by choosing a photon energy (and thus fcj_) 
that makes the ARPES spectral weight nearly zero for 
the upper AB sub-band and at the same time maximiz- 
ing the intensity of the lower sub-band. In principle then, 
the number of graphene planes that are part of AB pairs 
can be determined by counting the occurrence of these 
bi-layer bands relative to the number of linear bands from 
single graphene planes. 

While a straight forward procedure, the limitations of 
this method must be kept in mind. First, for typical pho- 
ton energies the photoelectron mean free path ranges be- 
tween l-5A).l^aH55lThis means that ARPES only probes 
the stacking order in the first 3-4 layers. Because the 
growth process of C-face films in the oven environment is 
not known, it is difficult to extrapolate surface values of 
Pab to the entire film. Also, merely counting cones and 
disregarding the emission intensity can be problematic. 
The intensity of a cone depends both on the depth of the 
graphene sheet in the stack and its relative area within 
the photon beam. Counting cones alone allows the possi- 
bility of overweighing small relative area graphene sheets 
at the expense of large but deeper sheets. This effect can 
be minimized by using as small a probe area as possible. 

With these provisos in mind, we have estimated Pab 
for a 9-10 layer graphene film using ARPES. To improve 
statistics, images similar to those in Fig. [T0| were collected 
on up to six different locations on 2 different samples. A 
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total of 110 cones were measured. Data was taken at 
a photon energy of 36eV to maximize the sensitivity to 
the lower bi-layer sub-band in Fig. [iBjb). We find that 
the number of planes in AB pairs is ^ 15% of the total 
planes in the film. Within a few percent, this fraction is 
the same for graphene planes rotated both (p = 30° and 
0° from the SiC (2130) direction, indicating that there 
is no preference for Bernal stacking for film rotated in 
either principle direction. While an estimate, the ARPES 
value of Pab corroborate the SXRD results and places the 
concentration of AB stacked planes to be less than 20% 
of the film. Both techniques support the conclusion that 
the rotated graphene sheets are interleaved in the film 
and that the rotational stacking is ordered rather than 
random. 



IV. CONCLUSION 

Epitaxial graphene grown directly on both the 
SiC (0001) Si- face and (0001) C-face has exceptional film 
quality.™^ These films are atomically flat with the 
graphene sheets being continuous over macroscopic dis- 
tances (if not the entire crystal surface). Furnace growth 
methods developed in the last 5 years are currently the 
best method to produce exceptionally large macroscopic 
graphene sheets. In fact STM studies have yet to find a 
single example of a discontinuous top layer in a C-face 
film. In this article, we show specifically how the unique 
structural and electronic properties of C-face epitaxial 
graphene continue to make this material a serious candi- 
date for graphene electronics. 

Unlike Si-face graphene where the stacking is AB or 
Bernal, '^■^ C-face graphene has an unusual stacking.l^ 
On this face the graphene sheets are not only oriented 
30° relative to the SiC (2130) direction (typical of Si- 
face graphene), but are also rotated with a distribu- 
tion of angles around 0°. The two orientations occur 
with equal likelihood. These rotated planes are inter- 
leaved in the film and are not part of isolated Bernal 



stacked films with multiple rotated domains like those in 
HOPG graphite. Furthermore, surface x-ray diffraction 
and angle resolved photoemission spectroscopy experi- 
ments show that the rotational plane stacking is highly 
ordered and not random. The exceptional large graphene 
sheet size and highly ordered rotational stacking of C-face 
epitaxial graphene demonstrate that this is a new form 
of graphene stacking and not the disordered soot known 
as turbostratic graphite.l^ 

The result of this type of ordered rotational stacking 
is that there are few adjacent planes in the film that are 
AB stacked. The concentration of these Bernal stacked 
planes is less than 20% meaning that they can be con- 
sidered as faults in an otherwise ordered rotationally 
stacked film. This unique stacking leads to a symmetry 
change in adjacent non-Bernal layers that makes these 
graphene sheets act electronically like a stack of isolated 
graphene sheet s.l^^lEol Because these non-Bernal stacked 
layers make up the majority of the film, nearly the en- 
tire graphene film beha ves elect ronically like a stack of 
isolated graphene sheet s.'^^'^^'^H 

While the reason for this type of rotational stacking 
is not understood in detail, the observed rotational dis- 
tribution shows a strong correlation with graphene-SiC 
substrate commensurate structures that indicates an im- 
portant orientational ordering mechanism at the SiC- 
graphene interface. 
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